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ABSTRACT 


We wish to determine how effectively a lubricant film prevents 
metallic contact between two rolling surfaces (such as in ball bear- 
ii^s) as a function of surface roughness parameters. The param- 
eters considered are the spectral moments of the two-dimensional 
surface obtained by superposition of the two rolling surfaces. We 
consider the peak height distribution, estimation of one -dimensional 
profile spectral moments, and the estimation of two-dimensional 
surface moments from several profile measurements. Also given is 
an asymptotic relation between the mean film thickness and contact 
occurrences. 
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TWO-DIMENSIONAL GAUSSIAN PROCESSES APPLIED TO 

THE DETERMINATION OF CONTACT BETWEEN 

LUBRICATED ROLLING SURFACES 

Steven M. Sidik 
Lewis Research Center 

W SUMMARY 

The basic problem under consideration is that of describing the 
behavior of two rolling surfaces separated by a thin film of lubricant. 
We discuss a procedure for determining the mean film thickness from 
experimentally obtained electrical contact measurements. We are 
specifically concerned with rolling surfaces such as occur in ball 
bearings. 

The ball surfaces are represented as Gaussian random processes. 
The development first considers the situation where contact occurs 
along a one -dimensional line segment. Both exact and approximate 
relations are developed between the mean film thickness and the 
amount of electrical contact. The development is then extended to 
the situation of contact over an area. An approximate relationship 
between mean film thickness and electrical contact is developed. 

The approximation presented involves the determination of the 
probability density function of peak heights on a random surface and 
the density of such peaks. 

One method of obtaining surface roughness information is to ob- 
serve relative surface height along one -dimensional profiles. We 
discuss the estimation of spectral moments of these one -dimensional 
processes. We then discuss a method of combining this information 
from several profiles in different directions to compute the spectral 
moments of the two-dimensional process. 
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We conclude with a brief discussion of the need for filterii^ the 
profile process in order to obtain better agreement with the contact 
process. 


INTRODUCTION 

The basic problem under consideration is that of describing the 
behavior of two rolling surfaces separated by a thin film of lubricant. 
This situation arises quite naturally in the study of ball bearings. 

The purpose of the lubricant is to separate the metallic surfaces in 
order to prevent bearing wear and thereby lengthen the effective life 
of the bearing. This report discusses a means of experimentally 
determining the mean film thickness from electrical contact measure- 
ments obtained by an experimental technique described in the first 
section. 

The proposed solution of this problem involves representing the 
ball surfaces as Gaussian random processes as was done in refer- 
ence 1. Supporting data for this is presented in Williamson (ref. 2). 
We begin the development by first considering the hypothetical situa- 
tion where contact takes place alor^ a one -dimensional line segment. 
That is, we consider two one -dimensional stationary and ergodic 
Gaussian random processes whose mean levels are separated by a 
lubricant film of mean thickness h. We then develop a relationship 
between h and electrical contact occurrences. An approximation 
is presented for the case when the mean film thickness is large com- 
pared to the surface roughness. This part of the development is 
similar to that of reference 1 except that an erroneous derivation in 
reference 1 is corrected. 

We then extend the study to the real problem where contact actu- 
ally occurs over an area. That is, we are now concerned with two 
two-dimensional stationary and ergodic random processes. This 
presents considerable difficulties in problem definition and analysis. 
An approximation is developed for the case when the mean film 
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thickness is large compared to the surface roughness. In order to 
achieve this result a number of results from and extensions of the sub- 
ject of "statistical geometry" as developed by Longuet -Higgins (refs. 

,3, 4, and 5) are presented. Also used are some results from Nayak 
(ref. 6), Cramer and Leadbetter (ref. 7), andKac and Slepian 
(ref. 8). 

The approximation developed provides an equation for the electrical 
no -contact time fraction as a function of the mean film thickness, Hert- 
zian contact area, and surface roughness parameters. This equation 
involves the probability density function of peak heights on a random 
surface and the number of peaks per unit surface area. 

One method of obtainii^ surface roughness information is to ob- 
serve relative surface heights along one -dimensional profiles. We dis- 
cuss three methods of estimating the spectral moments of the resulting 
one -dimensional processes. The electrical contact prediction is a 
function of the spectral moments of the two-dimensional surface ob- 
tained by superimposing the two rolling surfaces. Thus a method of 
combining information from several profiles in different directions to 
compute the spectral moments of the two-dimensional process is dis- 
cussed. 

We conclude with a brief discussion of the need for filtering the 
profile processes in order to obtain better agreement with the contact 
process. 

It is appropriate to remark at this point that we consider only a 
simplified version of the real problem. Some of the complications of 
real life are bearing spin, bearir^ sliding, viscosity effects, non- 
uniform lubricant thickness, and several others. A survey of elasto- 
hydrodynamic lubrication and related topics is presented by Tallian 
(ref. 9). 



4 


ONE -DIMENSIONAL ANALYSIS 
Relation of Film Thickness and Surface Roughness 
to Electrical Contact Occurrences 

The situation is represented in figure 1. There are two balls, con- 
veniently referred to as the upper ball and the lower ball. Under a con- 
dition of load, because of local elastic deformation of the ball surfaces , 
there is a region of contact where the surfaces can be considered as 
locally parallel. This region is called the region of Hertzian contact 
and is illustrated in further detail in figure 2. 

The authors of reference 1 have shown that, within the Hertzian 
area, the two surfaces can be closely approximated by two Gaussian 
random processes separated by an average distance h which repre- 
sents the average lubricant film thickness. The upper surface is de- 
noted Z^(x). We use the notation Z^^(x) ~ N(ju.^, C^(r)) to denote that 
Z^(x) is distributed as a Gaussian random process with mean and 
auto covariance function (acvf)C^(r). The lower surface is denoted by 
(x) and we have (x) ~ (r)). We assume both Z^ and 

Z^ to be stationary and ergodic. The distance between the surfaces 
at any point x is described by Z (x) - Z (x). Thus, we consider 

Vt ^ 

the composite process defined by 

{(x) = Z^(x) -Zj(x) -h (1) 

This process is illustrated in figure 3. We have (Bendat, ref. 10) 

^(x)~ N(0 ,C(t)) 

where C(t) = C^(t) + (t). Note that C(0) = the variance or 

mean square of the process. Whenever ^(x) < -h, the surfaces are 
in the state of metallic contact. Obviously, the physical process can- 
not have ^(x) < -h. We assume that |{x) < -h corresponds to a state 
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of elastic deformation of the surface asperities, and that the process 
of elastically deforming and reforming asperities has no appreciable 
effect on the characterization of the surfaces in contact by equa- 
tion (1). 

In order to determine, experimentally, when the surfaces are in 
metallic contact, it is possible to employ the electrical contact method 
described in reference 1. We describe this technique briefly as fol- 
lows. 

With reference to figure 1, we might apply a voltage across the 
upper and lower balls. Thus, whenever the balls are in metallic con- 
tact, a current flow will occur through the circuit. A continuous re- 
cording of the voltage will then provide a record of whenever the balls 
are in metallic contact. It is important to note that no voltage (i. e. , 
potential) will be recorded whenever there is metallic contact any- 
where within the region of Hertzian contact. A voltage will be ob- 
served only for those time intervals such that the lubricant film is un- 
broken throughout the region of Hertzian contact. 

We assume that the surfaces are rollii^ against each other at a 
constant velocity with no sliding or spinning. Then if a voltage - 
averaging device is used during a test, the time-averaged fraction 
of the applied voltage is a function of the distance -averaged fraction 
of rolled -over -distance which corresponds to the state of metallic con- 
tact. We denote the distance -averaged fraction of electrical contact 
as T.. The no-contact fraction is thus given by 1 - T . We wish to 
derive a relation between the average film thickness h and the con- 
tact fraction T^ so that experimental measurements of average volt- 
age may be used to estimate h. 

Derivation of Contact Fraction 


We begin by describing what happens during an electrical no- 
contact occurrence. Let the length of the Hertzian contact be denoted 
by d. For each excursion of ^(x) above the level |(x) = -h, electri- 



6 


cal contact is not broken until the left edge of the Hertzian contact in- 
cludes the first point at which ^(x) > -h. That is, electrical contact 
is not broken until the right edge of the Hertzian contact is d units of 
distance past the first point for which |(x) > -h. Electrical contact 
is reestablished as soon as the right edge of the Hertzian contact in- 
cludes the first point for which |(x) = -h again. Thus, if the excur- 
sion above -h is of length L, there is a loss of electrical contact 
only for the length (L - d)“^ where (x)'*’ = MAX { x, 0 } , 

Define the followir^ random variables (rv) associated with the 


random process 

l(x): 

«a(h) 

rv denoting ler^ths of excursions of | (x) above 
level -h 

N^(t,h) 

rv denoting number of crossings of | (x) from 
below -h to above -h in interval (0, t) 

N^(t,hl6sd) 

rv denoting number of crossii^s of |(x) from 
below -h to above -h in interval (0, t) where 
excursion is of length less than or equal to d 

N^(t,h|6>d) 

rv denoting number of crossii^s from below -h 
to above -h where the excursion is of length 
greater than d 

N(t,h) 

rv denotii^ total number of crossings of level 
in either direction 


Let the random variables 6^j(h), N^^(t, h), N|^(t ,h|6 ^ d), and 
Nb(t.h|e > d) be defined similarly for excursions and crossings from 
above -h to below -h. Cramer and Leadbetter (ref. 7, ch. 11) 
present rigorous proofs that these can, in fact, be treated as random 
variables. The following developments are presented in a heuristic 
manner but can be made rigorous through application of the results 
in reference 7. 
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We denote the distribution functions of 0^{h) and 6|j(h) by 
and Bj^( 6|^), respectively. Although there has been much research 
directed toward identifying these distributions, they are not known ex- 
cept for special values of h and certain forms of the acvf C(t) (e. g. , 
see Slepian, ref. 11). We assume the expectations of 6g^(h) and 5|jj(h) 
exist and denote them, respectively, as t^ and t^^. 

Of primary interest are the two conditional expectations and 
defined by 




= E 


= E 


^(h) 1 6^(h) > ^ 
«a0’) 1 5a(h) s ^ 


where d is the length of the Hertzian contact. Since 6g^(h) is a posi- 
tive rv, we have 


t.. = T^^Pr 
a N 


6a(h) 


> d 


+ TpPr 


5a(h) £ d 


For sufficiently large t, the following statements are true (~ de- 
notes approximately equal): 

2 

Na(t,h|6^(h)>d) ^^_ 

Ng^{t, h) 


6^(h) > d 


N_ t,h 6. 


^(h) s d) 


N^(t,h) 


Pr 6„(l>) s d 
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Because of the stationarity and ergodicity assumptions, the approxi- 
mate equalities approach equality as t — <». Likewise, the following 
derivations involving approximate equalities become exact in the 
limit as t - Thus, for sufficiently large t, the total no-contact 
distance is approximated by the number of excursions above the 
level -h which last longer than d, multiplied by the average amount 
by which such excursions exceed d. Thus, 

,. Na(t.h^(h)>d)(T^-d) 

t 

But we have 

Thus t 

^ Njt,h|6^(h)>d)(T^-d) 

(ta . y NM) 



( 2 ) 
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In order to utilize equation (2), we must know the distribution of 6g^(h) 
and the expectation of 6^(h). In general, these are unknown. There 
are some asymptotic results for 6j^(h) as h — <», Cramer and 
Leadbetter (ref. 7) also have a number of results concerning the expec- 
tation of the excursion ler^th which are exact. It must be noted that 
these exact results are, unfortunately, inapplicable since the quantity 
of interest is E(L - d)'*', not E(L). Equation (2) reduces to the ex- 
pression derived in reference 1 when the excursion lengths are assumed 
to be exponentially distributed. This assumption is open to question 
although it is true in the limit as h — >», 

Computational Method as h Approaches 

One of the primary purposes behind studying the EHD process is 
to design bearing systems where there is little or no metallic contact 
between surfaces. This implies that primary concern be given to the 
processes where h is large compared to the surface roughness. That 
is, where h/a is large. 

For this situation some results from Cramer and Leadbetter 
(ref, 7) and from Rice (ref. 12) and Kac and Slepian (ref. 8) are valu- 
able. These authors consider (1) the distribution of the lengths be- 
tween successive upcrossings of a h:^h level and (2) the distribution 
of the lengths of excursions above a high level. Since the Gaussian 
process is symmetric about the mean, it is evident that simular results 
apply to excursions and crossings below a low level. 

With regard to item (1), let F (x) denote the distribution function 
of the interval between successive downcrossings of the level ^(x) = -h, 
and let 1/0^ denote the mean of this distribution. It is shown in refer- 
ence 7 that 
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lim 
h — oo 



1 - 


e 


-X 


= Pr (the interval between one downcrossing and the next 
of |(x) = -h is ^/9^) 


Since 


e C"(0) ^-hV2 


we have 


Hence, the expected length between downcrossings increases rapidly 
and the stream of downcrossings becomes Poisson in nature. 

With regard to item (2), let F^(x) denote the distribution function 
of the lengths of excursions below the level ^ (x) = -h and let 6^ denote 
the mean of this distribution. Then it is shown in reference 7 that 

lim F (0 x) = 1 - exp 
h — 00 

= Pr (excursion length is 



Since for large h 


9 


e 


1 / 2 ^ 

hVc’COV 


we have 
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0 *-0 

h — oo 

Hence, the expected length of the excursion decreases rapidly. 

These two results may now be combined to give a limiting expres- 
sion for T„. In particular, for sufficiently large h, we will have ex- 
cursions below -h occurring quite infrequently; and each excursion 
causes an electrical contact whose duration will be, on the average, 
d+ 9 Thus, for sufficiently large h/a and t, 

Njt,h)(d + 0 ) 

m i-s-» ^ 


We note that Rice (ref, 13) has shown 


E 


Nb(t, h) 
t 



and hence 


Tc 



-C"(0) 


C(0) 


1/2 




Since C(0) = a , the variance of the Gaussian process, this equation 
may be rewritten and solved in terms of h/cr. In order to simplify this, 
we assume is small relative to d so that 





12 


or 






-C**(0) 

2 

_ O’ 
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TWO-DIMENSIONAL ANALYSIS 
Derivation of Contact Fraction 

At this point we make the observation that the Hertzian region is 
(for two balls) a circular area. The authors of reference 1 began 
their development with the assumption that the Hertzian region is an 
arbitrarily narrow and approximately rectangular area. They then 
developed what is essentially a one -dimensional theory. After con- 
sidering the one -dimensional theory, they then made a rather weak 
heuristic generalization to the two-dimensional theory. 

The results of Nayak (ref. 6) show that a naive analysis assuming 
that a profile trace in a single direction may be used directly is erro- 
neous. He shows how the two-dimensional moments may be estimated 
from a one -dimensional trace under the assumption of statistical iso- 
tropy. Both Nayak (ref. 6) and Longuet -Higgins (refs. 3, 4, and 5) 
consider such statistical descriptions as distribution of maxima and 
minima, distribution of peak heights, mean summit curvature, and so 
forth, Longuet -Higgins also discusses a method of determining a se- 
quence of estimating functions which converge to the two-dimensional 
spectrum. 

The statistical problem requiring solution before electrical contact 
analysis can be rigorously applied is now outlined. A complete solu- 
tion to this problem is yet to be developed. In fact, an appropriate 
formulation would be useful. 
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The composite process of interest can be described similarly to 
equation (1) as 


|(x,y) = Z^j(x,y) - (x,y) - h (4) 


where 


(x, y) are the coordinates of the reference plane 



Z ^ (x, y) 




and hence 


J(x,y)~ N^0,C(T^,Ty>) 


where 


C(VTy)=C„(T^,Ty) + C^(T^.Ty) 

2 

and C(0, 0) = a , the variance or mean square of the process. 

Figure 4 presents a visualization of the metallic -contact areas in 
the plane parallel to the (x, y) plane at level -h. For simplicity, we 
assume the x-axis corresponds to the direction of rolling. The area 
over which the Hertzian area passes as the surfaces roll over each 
other is bounded by the two horizontal dashed lines denoted y = 
and y = Y 2 . The areas for which ^(x, y) < -h are indicated by the 
enclosed dotted regions. For simplicity, we assume that the Hertzian 
area is circular. 

The left circle corresponds to the position of the Hertzian area 
when electrical contact is first broken. The right circle corresponds 
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to the position of the Hertzian area when electrical contact is first made 
after passing over the ”hill. ” 

From figure 4 we see that the quantity whose distribution we need 
to known is X. Since the distribution of X appears too complicated 
to derive at the current state of knowledge concernii^ two-dimensional 
random processes, we will consider a potentially useful approximation 
for the case when h/a is large. 

Computational Method as h Approaches «> 

As h/a gets large we might expect a similar asymptotic Poisson 
character cf the number of ’’valleys** to occur in the two-dimensional 
process as occurs in the one -dimensional process. This point should 
be verified. In this case, the metallic-contact occurrences are more 
appropriately illustrated by figure 5. This figure presents a visualiza- 
tion of the metallic -contact areas in the plane |(x,y) = -h. The areas 
for which 4(x,y) -h are indicated by the enclosed dotted regions. 

The area over which the Hertzian area passes as the surfaces roll over 
each other is bounded by the two horizontal dashed lines denoted y = 
and y - Y 2 . The dashed circle on the left of the region denoted I cor- 
responds to the position of the Hertzian area when electrical contact is 
first made. The point (Xj^j, y^) denotes the coordinates of the point 
where contact is first made. The dashed circle to the right of the 
region I corresponds to the position of the Hertzian area when elec- 
trical contact is finally broken. The point (xj^ , y^^ ) denotes the coor- 
dinates of the point where contact is first broken. We let X denote 
the difference X|^ - x^^^. If X is relatively small compared to the 
diameter of the Hertzian area, we may introduce the following approxi- 
mation. It is reasonable to expect y^ and y^ to be quite highly and 
positively correlated. Thus, assume y^ and y^ are perfectly corre- 
lated and that y^ is uniformly distributed over the vertical distance 
bounded by y = Yj^ and y ~ Y 2 . In this case we have 
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E(L) = 2E(L J + E(X) 

«2E(Lj^ 


K the shape of the Hertzian area and its dimensions are known, it is a 
simple matter to compute E(L|jj). 

Let D denote the expected number of summits per unit area 
sujtn 

and p(4*) denote the probability density function of the summit heights, 
I*, in the composite process defined by |(x,y). Since the area over 
which the Hertzian region passes is of width w = - Y 2 , we obtain 

the expected number of summits occurrir^ within the strip per unit 
distance along the x-axis as The probability that two peaks 

occur within the Hertzian region is negligible by the Poisson nature of 
the peak occurrences. The average fraction of contact distance T^ 
is thus given by multiplying the expected number of summits that ex- 
ceed the level h per unit distance along x, times the average distance 
contact is maintained for each excursion. That is 



(5) 


Peak Height Probability Density Function 

The density function p(^*) is really a conditional density of the 
surface height conditioned on the event that the point is a summit. 
Since the conditioning event has zero probability, great care must be 
exercised in the definition of the conditional probability density. This 
has been discussed by Kac and Slepian (ref. 8). 

Nayak (ref. 6) presents a derivation of p(^ ). The derivation 
presented here shows that the results correspond to a ’’horizontal 
window” conditional density as defined by Kac and Slepian. This is 
quite important, as only ’’horizontal window” arguments lead to the 
pdf possessing the ergodic property. Let 



16 


^4 = 

«i = — l(x,y) 

3x^ 

^2 = -®— |(x,y) 

ax3y 

^3 = -^^(x,y) 
0 2 
3y^ 


( 6 ) 


^5 = — ^(x,y) 

^ 0X 


^6 = -^(^’y) 

^ 3y 


Since ^ (x, y) is assumed to be a stationary and ergodic Gaussian 
process, there is an autocorrelation function R(t , t„) and a spectral 

n ^ y 

density function 6? {<p^, cp^) associated with the process. The function 

R(t„, T.J is defined by 
A y 

V " = -i E + T^, ty + Ty){(t^, ty)j 

A duality between R(Tj^, t ) and e/ <p ) is defined by 
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R(t^, 





c/ {<p^, <p^)d<p^A<p^ 


and 


e/ {9^, <py) = 


^00 ^00 


-00 -00 


-i{r (p -VT (p ) 

e y y R(Tjj, Ty)dT^dTy 


The spectral moments of fi/ are defined by 



(p\<p\ CS {<py.y (p^)d(p^d(p^ 


a) 


In terms of these spectral moments, it is well known that 4 ~ N(0, S) 
where 
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To derive the distribution of peak heights, consider the area of the 
x-y plane dA = (x,x + 6^), (y, y + 6y) where we have at some point 
(Xo> Yq) in dA 


«3 = «3(Xo,yo) = 0 

l2 = ?2<’‘0>i'0> 

4i{3-{2^0 

The region where the restrictions on ?3 are satisfied 

is denoted V. 

We wish to define the conditional density function 

P(^ 4 ll 5 = = summit) 

We proceed with a ’’horizontal window” limiting argument as presented 
in Kac and Slepian (ref. 8) since we desire ergodic properties of this 
pdf to hold in the plane (x,y). Thus 

p(5^[ summit) 


5 5 " ^ 5^*0’ ® 

«6=?6<*0>5'o)= ° 


lim 

dA - 0 


P P Pr ^Ug ^ 

/ / / \ / . / P(^i.-4e)^5^6P^P^2^3 

1^5 ^6 J 
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where Ug,Ug), (ig, Ug) are lower and upper limits on and |g im- 
posed by the conditions 

I g « 0 in dA, $ g(Xg, yg) = 0 

and 

4g « 0 in dA, e6<^0’ V = ® 

Now ^g=3|/3x and B^/dy so that 

where 

|J| = Ul«3 -«2l 

Thus p(4^ I summit) is also given by 



The limit is undefined as 6 , 6. — 0 since the ratio approaches 0/0, 

^ y 

Use L’ Hospital’s rule with respect to 6^ and 6y, noting that 
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ae^asy 



0 '^0 


p(^i- . 46)^^y 


-•(0, 0) 


lim 



0 pU X’ 5 2’ ^3’ ^4’ ^x’ ^y^ ” X’^2’^3’^4’^’®^ 

0 


to obtain 


p(^ -[ summit) = 


jH J 1 p(e 1> ? 2- « 3’ 4 4’ 1^ 2‘*« 3 



l•I|p({l.l2.{3.l4.0,0)d^Jd{2d^3 





« 1?3 -«2 P (« 1 - 52 ’« 3 '« 4 ’®’“)^ 1 ^ 2‘^3 


Note that 


P(4 x»? 2’^3’^4’ 


1 -Q/2 

— - — e 

^^2^4 


where 
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Q - Ui>42»?3»?4) 




(10) 


^2 " ”' 02*”20 



= det S 




-i 

I 

“oo 


Nayak (ref. 6) has also presented an argument which concludes that 



And hence, the normalizing factor in equation (9) is simply 
The only problems remaining are the evaluation of and 



To this end, let 
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Usii^ some results about inverses of partitioned matrices (Graybill, 
ref. 14) 



Then the quadratic form Q of equation (10) is given by 
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■" ^ 8 ^ 12^4 ^ 4 ^ 21^8 ■*■ ^ 4 ^ 22^4 


= ^ s^ii^s + ■*■ ^4^22^4 


= 4 gRll 4 s ■*■ 2 ^ 8^11 ^11^12^4) (^4^21^11) 

X (RiiI^12^) + ^^22^4 

■ (^4^21^11)^11(^11^12^4) 

= (^8 + ^11^12^4) ^11 ^8 + ^11^12^) 

+ ^4 (^22 ■ ^21^11^12)^4 

(let z= 


since 


= (Is - 2 ) F - 


ll 


m 


00 


(■ 


,-1t 


^22 " ^21^11^12 





“^00 
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where c 
And 

-R 


and hence 
Thus 

For fixed 
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11^12 ■ 


1 T 
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00 


-F-lj 


= - [f . (cF-^O 
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00 


= - fci 
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00 


= + — 
^00 


^4 


I -z=l - —^l 
00 


m 




1 * 

m 


Is Z ) P'^ - cF'^Zz'^F"^) Z 


00 


1 * 

m 


00 


(14) 


we have 


N(0,F) 
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and hence 


”^00 \^00 / 

The determination of the pdf of the height given we are at a summit 
thus reduces to the problem of evaluating 




e d| 


(15) 


where 


y- = 



*”oo 


'^20 

-mu 

-”^02 


(16) 


There is, in general, no closed form expression for evaluating this 
integral. This can be seen from equation (15) which shows the integral 
to be related to the expectation of a function of normal r. v. *s over a 
conical region of the space. In fact, the particular conical region is a 
cone generated by a ray rotated around the line ? 4 3 with ^ - 

and Ig both negative. The angle of the ray with 4 = ^3 is 45®. 

From this representation we can see that the point i± will intro- 
duce asymmetry. There is a special case of the process we are con- 
sidering which retains enough symmetry so that the integral can be 
evaluated as a closed form. The special case is when the surface is 
isotropic. 
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Isotropic Processes 

An isotropic random surface is one in which the probabalistic be- 
havior of the surface, alot^ a profile observed at an angle 9 , is inde- 
pendent of 6 . In particular, this implies that 





27 


Under the isotropy assumption, the integral of equation (15) can be 
evaluated analytically (e.g. , Nayak, ref. 6) to obtain 
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It may be noted that the properties of an isotropic random surface 
may be estimated from the statistics of a single profile measured 
alor^ any direction in the plane. This is not true for non-isotropic 
surfaces, however. 

The surface of a ball may reasonably be assumed to be isotropic 
because of the way in which they are finished. Hence, when two balls 
are in contact, the composite process will also be isotropic. This is 
the case considered in reference 1 and Nayak^s (ref. 6) results are 
applicable. 

From the way in which ball bearing raceways are finished, how- 
ever, one would expect the surface to have grooves and ridges 
elongated in the direction of grindii^. This can be visually verified 
in some instances. The composite surface process resulting from a 
ball rolling in a raceway is unlikely to be isotropic. To allow for analy- 
ses involving balls rolling in raceways, we now consider methods for 
estimating the spectral moments m^^ so that equation (5) may be used. 

DETERMINATION OF SPECTRAL MOMENTS 

The use of equation (5) to estimate T depends upon the ability to 

jjc ^ 

numerically evaluate a>nd p(| ). These quantities, in turn, de- 

pend upon the spectral moments m^j of the composite process defined 
by the sum of the individual surfaces. We outline a method of samplii^ 
the individual surfaces of the ball and raceway aloi^ selected prcrfiles 
obtained by use of a profilometer. We begin by discussing methods of 
estimating the moments along one -dimensional profiles. We then use 
certain of the results of Longuet -Higgins (ref. 3) which relate profile 
one -dimensional moments to the two-dimensional moments of the en- 
tire surface. 


Estimation of Profile Moments 

As indicated, the actual surfaces are sampled in a one -dimensional 
manner with a profilometer. The usual method is to draw a fine- 
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pointed stylus across the surface and obtain an analog record of the 
relative surface height along the profile. This analog record may 
then be digitized and analyzed on a digital computer. There are three 
possible methods of analysis. 

1. Numerically compute the sample one -dimensional spectrum by 
use of a Fourier transform method (refs. 15 and 16). The moments 
of the spectrum may then be numerically calculated by quadrature. 
This method is relatively straightforward and requires no further 
discussion. 

2. Count the number of relative extrema (maxima and minima) 
and the number of zero -crossings of a profile sample. K we let 

D_ = density of zero -crossings 

Dg = density of extrema 
2 

a = estimated variance of process 

Then it can be shown that estimates for the moments can be defined as 


and 


2-2rx2 
”2 = ” " 





(19) 


Naynk (ref. 6) presents a more complete discussion of this method and 

^2 

also indicates a method by which D^, and cr can be obtained 
directly from the analog relative surface height record by electrical 
devices. 

3, Assume some physically realistic but tractably simple family 
of analytical models for the autocorrelation function. Numerically 
estimate the autocorrelation function from the digitized profile record 
via Fourier transform methods. Find the best -fitting member of the 
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family of autocorrelation functions (i. e. , estimate the unknown param- 
eters) and estimate the moments via the well-known relationships 
(Cramer and Leadbetter, ref. 7) 


m 



m„ = ^R^^)(0) (20) 

^ -2 
(J 

m4 = -lR<4)(0) 


This method requires further development. 

Much of the book by Cramer and Leadbetter is devoted to the con- 
ditions under which these moments (and, hence, derivatives of R(t) 
at T = 0) exist, and the relation of these moments to analytical prop- 
erties of random processes. Of particular interest is the relation of 
the moments to such things as the density of zero -crossings and the 
density of extrema. We refer the reader to Cramer and Leadbetter 
for full discussion of these topics and limit ourselves to one principal 
result. This result is that the i'"'^ spectral moment exists if and only 
if R^^\t) exists and is finite at t = 0. If the second spectral moment 
is infinite, then the expected number of zero crossings in any finite in- 
terval is infinite. 

It is thus interesting to note that some authors (refs. 17 and 18) 
have proposed exponential or exponential -cosine autocorrelation 
functions as realistic models for surface roughness. Bendat (ref. 10) 
has also proposed these as beii^ realistic in many other applications. 
The major theoretically desirable property of a Gaussian random pro- 
cess is that an exponential autocorrelation function implies that the 
process is Markovian. 
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Consider the following family of autocorrelation functions 

R(x) = e'l^l°' (21) 

where of = 1 yields an exponential correlation function. 

Then 



The various values of R^^^(O) and R^^\o) are given in the following 
table. From these values, we see that for or > 2 the expejcted number of 
crossii^s at any finite level in any finite interval of time is zero. If 
Q! < 2, then the expected number of crossings of any finite level in any 
finite interval of time is infinite (or undefined, but undefined in such a 
manner that we would expect the number of crossings to be infinite). 

These properties are difficult to accept from physical considerations. 


OL 

R^^^(0") 


R^^^(O) 

a>2 

0 

0 

0 

a= 2 

0 

0 

-2 

2 > a > 1 

0 

0 

-OO 

a = 1 

+1 

-1 

Undefined 

or < 1 

+00 

-OO 

Undefined 
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In order to resolve this point we propose the following class of 
models. Let 


p(t) = 1 + ajT^ + a2T^ + . . . + (24) 

q(r) = 1 + + bgT^ + . . . + b^r^”^ (25) 

where 

m - n = k > 0, b. > 0, b^^ > 0 

and 

R(r) = (26) 

q(r) 

Doob (ref. 19) has considered a similar approach to the definition of 
spectral densities. This representation seems reasonable as an ap- 
proximation to the exponential -cosine class for the following reasons 

1. Since p(r) and q(r) are both even, we have the required sym- 
metry. 

2. Since the degree of q is larger than the degree of p, and q 
is always positive, R(t) will be damped and R(t) — 0 as t 

3. If n is sufficiently large, and the roots of p are equally 
spaced, R(t) will have the same qualitative behavior as an exponen- 
tial cosine function. 

4. Since R^^^(O) = 0 for j ^ 2m, all spectral moments of order 
2m or greater are equal to zero. 


Multiple Profile Analysis 

Once the spectral moments of the profiles of the composite surface 
are obtained aloi^ several profiles we wish to use them to estimate the 
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two-dimensional moments. Xet 

m o = the n^^ spectral moment of the one-dimensional 
n, p 

profile taken at an angle of 9 from the x-axis of 
the (x, y) plane. 

Longuet -Higgins (ref. 3) has shown the following relationship to be true 

n 

mn_i ^cos*'‘V)sin^(e) (27) 

i=n 

Suppose that n profiles are taken at the (distinct) aisles 0., 
i = l,n. From these profiles we thus obtain 

niQ ~ ^00 ~ (28. 1) 

“ “^20^il ^ll^i2 ”^02^i3 

0. ^ ”^40^il ”^31^i2 ”^22^i3 ”^13^i4 “^04^i5 ^ 

(28. 3) 

Nayak (ref. 6) claims that three profiles taken at three distinct directions 
is sufficient information to determine the quantities and p(^ ). 

This seems unlikely to be true for the following reason. 

Longuet-Higgins (ref. 5) has shown that many of the gometrical prop- 
erties of a random surface are determined by the so-called invariants. 

Zero^^ order invariant: 



^0 " ”^00 
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Second order invariants: 

Sj = mg2 + ^20 

S2= “20“02 -“ii = det|s| 

Fourth order invariants: 

^ ”^40 ^*”22 “04 

Pj - >»04“40 ■ ^“13”31 + ®“22 

^3 “ <“40 + “22)<“04 “22> ' <“31 + “l3>^ 

= det F 

Since there are four invariants of order four, it may be possible that 
Dgum ^ expressed in terms of these four invariants. 

In such a case it would be sufficient to take four profiles and solve equa- 
tion (28.3) for the two-dimensional moments or the required invariants. 

It is clear, however, that five distinct profiles will be sufficient since 
there would then be five equations to be solved for the five unknown fourth 
order two-dimensional moments. The five equations for the zero*^" and 
second order moments may be solved by least squares. If more than five 
profiles are obtained, then the fourth order moment equations may also be 
solved for by least squares. 

Prefiltering of Signals 

Whitehouse and Archard (ref. 17) have discussed briefly the presence 
of very low and very high frequency components of the spectrum and the 
resultant effect upon surface analyses. This deserves some further dis- 
cussion in relation to the theory of electrical contact occurrences as de- 
veloped herein. 
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We first consider high frequencies. A typical profile that might result 
from a surface which contains primarily low-frequency components super- 
imposed upon a signal of high frequency with small amplitude is illustrated 
in figure 6. It may easily be recognized that the high-frequency component 
contributes very heavily toward the number of ’’summits” contained in the 
signal. Not all these ’’summits” have physical s^nificance with relation 
to electrical contacts, however. We may assume that since the asperities 
on the surface undergo eleastic deformation, the small asperities will dis- 
appear when the larger asperities come into contact. Thus, the data re- 
corded for description of the surface should be filtered so as to remove 
the high-frequency component. 

Figure 7 presents a typical profile that might result from a surface 
which contains primarily medium -frequency components superimposed on 
a signal with low frequency and low amplitude. In this case the low- 
frequency component will be reflected as very long and smooth ’’hills" 
and ’’valleys. ” But if these are quite large with respect to the Hertzian 
area, it is more reasonable to assume that the bearings will roll up over 
the hills and down through the valleys rather than cause very long contact 
and no -contact occurrences. Thus, the data recorded for description of 
the surface should also be filtered to remove such low-frequency compo- 
nents. 

Since there are little data available in sufficient detail, it is not clear 
that this problem need arise in practice. The most serious problem is 
with the high-frequency signals, since they contribute most to the number 
of summits. Whitehouse and Ar chard {ref. 17) present some evidence that 
at least some surfaces may be described as in figure 6. In addition, they 
also comment that in some instances the smaller peaks undergo plastic 
rather than elastic deformation and hence disappear rapidly with use. 

This lends support to the procedure of filtering out the high-frequency 
components and also presents a potential means for determining where 
the high-frequency cutoff point should be. 

In any event, what may be considered ’’high’’ frequencies and ’’low’’ 
frequencies is a subject for further research. 
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CONCLUDING REMARKS 

The basic problem under consideration is that of describing the be- 
havior of two rolling surfaces separated by a thin film of lubricant. We 
discuss a means of determining the mean film thickness from experi- 
mentally obtained electrical contact measurements. We are specifically 
concerned with rolling surfaces such as occur in ball bearings. 

The ball surfaces are represented as Gaussian random processes. 

The development first considers the situation where contact occurs along 
a one -dimensional line segment. Both exact and approximate relations 
are developed between the mean film thickness and the amount of elec- 
trical contact. The development is then extended to the situation of 
contact over an area. An approximate relationship between mean film 
thickness and electrical contact is developed. 

The approximation presented involves the determination of the proba- 
bility density function of peak heights on a random surface and the den- 
sity of such peaks. 

One method of obtaining surface roughness information is to observe 
relative surface height along one -dimensional profiles. We discuss the 
estimation of spectral moments of these one -dimensional processes. We 
then discuss a method of combinii^ this information from several profiles 
in different directions to compute the spectral moments of the two- 
dimensional process. 

We conclude with a brief discussion of the need for filtering the pro- 
file processes in order to obtain better agreement with the contact process. 
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Figure 1. - Idealization of lubricated balis in contact. figure 3. - Composite prxess defined by ^(x) = z^lx) - z^(x) - h. 

(Scale of Hertzian contact greatly exaggerated. I 



Figure 4. - Typical break -contact occurrence. 



□ No metallic contact 
Metallic contact 



Figure 5 , - Typical make-contact occurrence when contact areas are small. 



Figure 6. - Typical profile with high frequencies of small amplitude superim- 
posed on a (ow-frequency signal. 



figure 7. - Typical profile containing primarily medium -frequency components superimposed on a 
signal with low frequency. 


